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Abstract 

We investigate a spherical overdensity model for the non-clustering dark energy (DE) with the 
constant equation of state, (Jde in a flat universe. In this case, the exact solution for the evolution 
of the scale factor is obtained for general oj^e- We also obtain the exact (when uj^e = — §) and the 
approximate (when ui^e ^ ~\) solutions for the ratio of the overdensity radius to its value at the 
turnaround epoch (y = -^-) for general cosmological parameters. Also the exact and approximate 
solutions of the overdensity at the turnaround epoch £ = Pcl " stcr are obtained for general tu^e ■ 
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Thus, we are able to obtain the non- linear overdensity A = 1 + <5nl = C(|) at any epoch for 

the given DE model. The non-linear overdensity at the virial epoch A v i r = cTt^I is obtained 
by using the virial theorem and the energy conservation. The non- linear overdensity of every DE 
model converges to that of the Einstein de Sitter universe A^ S = 187r 2 + |) — 147 when z v - ir 
increases. We find that the observed quantities at high redshifts are insensitive to the different 
ujde models. The low-redshift cluster (z w - lv ~ 0.04, i.e., z ta ~ 0.7) shows the most model dependent 
feature and it should be a suitable object for testing DE models. Also as 0^ increases, the model 
dependence of the observed quantities decreases. The error in the approximate solutions is at 
most 2 % for a wide range of the parameter space. Even though the analytic forms of y and £ 
are obtained for the constant lj^, they can be generalized to the slowly varying u;d e . Thus, these 
analytic forms of the scale factor, y, and £ provide a very accurate and useful tool for measuring 
the properties of DE. 



1 Introduction 



Galaxy clusters are the largest virialized structure formed by the gravitational collapse of small 
density perturbations. The spherical collapse model (SCM) is a simple and powerful tool for 
investigating the evolution of nonlinear structure and bound systems in the Universe [1, 2]. 

Since the formation of the large scale dark matter potential well of clusters depends only on 
gravitation, the property of dark energy (DE) and the cosmological parameters can be constrained 
from the growth of large scale structure and the abundances of rich clusters of galaxies. There 
have been numerous papers investigating this problem with the spherical collapse model [3, 4, 5, 
6, 7, 8, 9, 10, 11, 12]. 

In order to prohibit the spherical overdensity from collapsing to a singular point in a finite 
time, we need to use the virialization in the SCM. The common mistake in literatures is that 
one uses the correct value of the finite virial radius of the overdensity from the virial theorem 
but still uses the incorrect value of the background scale factor for an Einstein de Sitter (EdS) 
universe. However, this leads to an incorrect conclusion that the bound systems are formed when 
the nonlinear overdensity is about 178. Equally, this leads to the incorrect critical linear threshold 
density 1.68 [13]. When we use the correct values of the virial radius and the scale factor at the 
virial epoch, these values are lowered to 147 and 1.58, respectively. 

In the next section, we review and analyze the SCM model in a flat universe including the 
non-clustering DE with the constant equation of state aide- We obtain the exact solution for the 
evolution of the scale factor (a) in terms of the ratio of it to its value at the turnaround epoch 
(x = ^). We also obtain the exact (when oj de = — §) and the approximate (when u> dc ^ ~\) 
solutions for the ratio of the overdensity radius to its value at the turnaround epoch (y = -^-) 
for general cosmological parameters. Also the exact and approximate solutions of the overdensity 
at the turnaround epoch £ = pc ^ stor are obtained for general Wdc- The given approximate 



= 2ta 



solutions have the less than 1.5 % errors over the parameter space in the concordance model. In 
Sec. 3, we use the virial theorem to obtain the correct values of x and y at the virial epoch. We 
compare the values of the virial epoch and the nonlinear overdensity for the different models. We 
conclude in Sec. 4. 



2 Spherical Collapse Model 



Background evolution equations of the physical quantities in a flat FRW universe with the matter 
and the DE are given by 
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Pm + 3Qp m = 0, (2.3) 

/ 9 de + 3(l+to dc )(^)p dc = 0, (2.4) 

where a is the scale factor, u; dc is the equation of state (eos) of the dark energy, p cr is the critical 
energy density, and p m and p dc are the energy densities of the matter and the dark energy, 
respectively. We consider the constant u; dc models only. 

We investigate a spherical perturbation in the matter density [1, 2]. /^luster is the matter 
density within the spherical overdensity radius R. The flatness condition is not held because of 
the perturbation in the matter. Thus, we have another set of equations governing the dynamics 
of the spherical perturbation [3] : 
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where r = 3(1 + u de ) (- - -) p dcc with < a < 1 , (2.7) 

where p dec is the DE density within the cluster, a = gives the clustering behavior of DE and 
a = 1 means the homogenous DE satisfying Eq. (2.4). For a = 1 case, the system does not 
conserve energy as it collapses from the turnaround to the virialized state. We can solve Eq. (2.7) 
for the constant a, 

I \ -3(l+w do )a / \ -3(l+w de )(l-a) 

p dcc (a,R)=p£ c — — , (2.8) 



where a ta and R ta define the scale factor and the radius at the turnaround redshift z ta - There is 
an ambiguity of p dcc in Eq. (2.8) because it is a function of both a and R except when a = 1 or 
0. Independent of a, the radius of the overdensity R evolves slower than the scale factor a and 
reaches its maximum size i? ta at z t & and then the system begins to collapse. 

Now we adopt the notations in Ref. [5] to investigate the evolutions of a and R, 

x = — , (2.9) 
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y = jt- ( 2 - 10 ) 
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Then from Eqs. (2.1) and (2.5), we obtain 
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where dr = H ta ^/Q m {x ta )dt, ( = ^^U a , Qta = 



(1 + Zta )3-3(l+^ dc )a 



\RtJ 



(5m 

-3(l+u, de )(l-a) 



Pm 



;i + z ta 



-3uj de 



Qcta — 

Hdcc 



, and x ta = 1 from Eq. (2.9). and 0[| e represent 



the present values of the energy density contrasts of the matter and the DE, respectively. Q ctSL 
becomes Qta when a = 1 (non-clustering DE). Also when uj^c = — 1/3 or —1, Eq. (2.12) becomes 
independent of a. We can solve Eq. (2.12) analytically in these cases. We limit our consideration 
to only the non-clustering DE. 

We are able to obtain the analytic solution of Eq. (2.11) 
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where F is the hypergeometric function and we use the boundary condition x = when r = 

(see Appendix for details). This solution is the same as the one in Ref. [13] when w de = -I. The 

above equation gives the well known relation a(t) oc sinli3 |^/j7[j c Hot when uj^c = — 1- From this 

equation we are able to find the exact turnaround time r ta for general Wde models, 
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where we use the fact that Xt a = 1- It is useful to have the exact analytic form of the turnaround 
time to investigate the accuracies of both numerical and analytical calculations. 

As expected, r ta depends on u& c , f^ n , and z ta as given in Eq. (2.14). We show these properties 
of T ta in Fig. 1. In the left panel of Fig. 1, we show r ta as a dependence on the values of O^, for 
the different DE models. We choose z ta = 1 in this figure. The dot-dashed, solid, dotted, and 
dashed lines (from top to bottom) correspond to oj^e = —1.2,-1.0,-0.8, and —1/3, respectively. 
If we have the same $7° n for the different models, then it takes longer time to the turnaround 
for the smaller values of code- The smaller value of Wde means the higher pressure against the 
gravity. Thus, it takes longer time to reach to the maximum radius of the overdensity for the 
higher pressure universe when we have the same amount of Q^. We also need more matter to 
overcome the higher pressure. We also show the Zt & dependence of Tt a for the different DE models 
when we choose Jl^ = 0.3 in the right panel of Fig. 1. Again the smaller values of lu^ c give the 
larger r ta . An interesting feature in this figure is that as zta increases so does Q t &- For large z ta , 



Qta 1 approaches to and thus r ta approaches to |F 
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Now we work on the evolution of y. We only consider the a = 1 case and Eq. (2.12) becomes 

1 



dT 2 



Cy- 



+ (l+3u Jdc )—x-^ 1+ ^y 

<Vta 



(2.15) 



For general u>de, there is no analytic solution of Eq. (2.15). However, we are able to obtain the 
exact analytic solutions for u>de = —1/3 and —1. First, we consider the case when u>de = — I- 
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Figure 1: r ta for the different values of uj^ e . a) r ta v.s. for the different values of u>d e = 
-1.2,-1.0,-0.8, and —1/3 (from top to bottom) when zt& = 1.0. b) Tt a v.s. z tSu for the same 



models when we choose = 0.3. 



In this case, the second term in the right hand side of Eq. (2.15) disappears and the equation 
becomes [14] 
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The solutions of the above equation are given by (see Ref. [13]) 



C 2 ± T . 



(2.16) 



ArcSin[^/y] - ^y(l - y) = VC13T ,when r < r t 

\Jy{ 1 -y) 



ArcSin[^/y] + - 



Ci3(t - T ta i 3 ) , whenr > r ta i3 , 



(2.17) 
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where £13 and r ta i3 denote respectively £ and r ta for uo& c = — ,-. Even though the above solutions 
given in Eqs. (2.17) and (2.18) are only true for uj& e = — w, these solutions can be good approximate 
solutions for general when we consider the high-redshift clusters. It is because the high redshift 
clusters which are formed at the high z v ; r require the high z ta in each model. As Zt& increases, 
so does Q t3i for the fixed value of ~ 0.3. Thus, the second term in Eq. (2.15) is negligible 
compared to the first one and this equation approximately equals to Eq. (2.16). Also these 
solutions are the same as the ones in the Einstein de Sitter (EdS) universe (f^ = 1 and f2[j = 0). 
It is because in the EdS universe, Q^ 1 = and the evolution equation of y is the same as in 
Eq. (2.16). Also the solution given in Eq. (2.17) is a good approximate solution for general lv^ c 
when we consider the evolution of y at early epoch. As y — > 0, the second term in Eq. (2.16) is 
negligible again. 



As long as we know the exact virialized epoch z v 
( for the different values of and z ta when u)& e = 



Eq. (2.18) provides the exact value of 
. We will investigate z v ir based on the 
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x y 
Figure 2: Both numerical and analytic solutions, denoted by dots and solid lines respectively, 
a) Cosmological evolutions of r for uj& c = — ^ as a function of x when we choose zt a = 0.8 and 
= 0.3. b) Evolutions of T ELS cl function of y for the same values of zta, and as in a) . 



virial theorem later. However, we first consider the simpler estimation to obtain £ values for 
the general u^ c . If we define z c as the redshift at which the cluster formally collapses to R = 
according to a spherical solution, then this collapse time is twice the turnaround time t c = 2i ta 
(i.e. t(z c ) = 2r(z ta )) [5]. If we adopt this fact, then we are able to obtain ( from the above 
analytic solution in Eq. (2.18) 
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well-known value of ( for the EdS universe [2, 15]. The general value of £ for lu^c = — 5 when 
7^ 1 is given by Eq. (2.19). After we obtain the value of £13 given in Eq. (2.19), we are able 
to find the values of x and y at any T ta < r < t c from Eqs. (2.13) and (2.18). We constrain 
t > T t a to show the evolution of x after the overdensity radius reaches its maximum. In the left 
panel of Fig. 2, we show the evolution of r as a function of x when zt & = 0.8 and fi^ = 0.3. The 
numerical solution (denoted by dots) is definitely the same as the analytic solution (solid line). 
In the right panel of Fig. 2, we show the evolutions of both analytic and numerical solutions of 
y(r) for the same values of z ta and as in the left panel of Fig. 2. 

Now we consider Eq. (2.15) when to^e = — 1- The differential equation of y for u^e = — 1 
becomes 
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We obtain c x = ( + Q^ 1 from Eq. (2.20) by using the initial condition ^| Zta = 0. The solution of 
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Figure 3: a) Cws>Ci3> an d Csk (from top to bottom) versus $7^ when we choose z ta , = 1. b) (s 
versus z ta when we choose 0^ = 0.3. 

the above equation is given by (see Appendix for details) 

c F ^EllipticF[ArcSin[B],c] + c P EllipticPi [A, ArcSin[B], c]^ = c 2 - r ,when r > r ta , (2.21) 

where EllipticF and EllipticPi correspond to the elliptic integral of the first kind and the incom- 
pelete elliptic integral of the third kind, respectively [16]. Also cp = 1 2 V 2 Qt a = an( j 



cp = 2 (QtaC 2 ) fT; ne integral constant c? is able to be obtained when we use the boundary 

(3 + Vl+4QtaC) 

condition y(rt a ) = 1. Even though this is the exact analytic solution of y for lij^o = — 1, this 
formula has some problems. First, this is inconvenient to use. Second, this formula does not 
produce the correct value when B > 1. However, it is useful to have this analytic solution when 
we derive the form of fitting formulae for £ and y for general Ude cases. 



C13 given in Eq. (2.19) is exact only for lo^ c = — g. Thus, we need to find £ for general uide- As 
we show in both uj^e = — § an d —1 cases, ( is a function of both u>^e an d Qta- Also it depends on 
a. Thus, we are able to derive a reasonable fitting formula for £ as a function of cu^e, Qta, an d a. 
When we investigate the numerical solutions of Eqs. (2.12) for general uid e , we obtain the fitting 
formula of £ 

£ k _ ^^^ 2 ^-0.724+0.157C m ta+a(l+Wdo)(l+3w de )(0.064-0.368n mt a) ^ ^ 

where r2 m t a = ^m(^ta) = ; — Tpr a nd it becomes 1 when £7^ = 1. Thus, this choice of the fitting 
formula is consistent with the fact that £ = (^-) 2 in the EdS universe. Also f2 m t a is a function of 
w dej z ta, an d f2^, and thus it is a function of all the variables in the consideration. Even though 
we consider a = 1 model in this paper, the above fitting formula is good for any value of a. The 
special choice of this fitting form of £ is independent of a when uj^e = — 4 or —1 in order to be 
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consistent with real models. In Ref. [5], C is given by another form for the fitting formula, 



'0.79+0.26n mta -0.06tJ de 
mta 



(2.23) 



In Fig. 3, we compare the values of Csk and Cws for uj Ac = — | with the exact C13 given in Eq.(2.19) 
for the different values of and z ta . We use z ta = 1 and normalize (s with (^) 2 in the left 
panel of Fig. 3. The dashed line depicts Cws and the solid and dot-dashed lines describe C13 and 
Csk) respectively. £i3 and Csk are too close to each other and look like a single line. We show £s 
as a dependence on z ta when we choose = 0.3 in the right panel of Fig. 3. Again both £13 
and Csk are very close to each other in most of the z ta range. We will show later that the fitting 
formula given in Eq. (2.22) shows the better fitting than that of Eq (2.23) for general Wde- 

As given in Eq. (2.15), y should be a function of cude, Qta 5 and Even though the function 
of y for the general value of tu^c should not be the same as Eq. (2.18), we modify this equation to 
obtain the approximate analytic solution of y as 

^ v 2.7+0.in mta -3.4(l+u; d e)-0.02Csk+3(zt a -0.6) 

~^mta I 
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This approximate solution is obtained from interpolating y values for the wide ranges of the cos- 
mological parameters and one can further work on the similar fitting form by changing interested 
ranges of the cosmological parameters. There are several reasons for the specific choice of the cor- 

/ \ 2.7+0.in mt a-3.4(l+w de )-0.02< s k+3(>ta-0.6) 

rection term I — ( 1+ ^ dc ) Q mta ) in the above equation. First, 

we know that the first two terms in Eq. (2.24) are exact when ojd c = — j- Thus, the correction 
term should not appear for lo^c = — | and thus we put (1 + 3o;de) in this correction. Also as zt a 
increases, the contribution of the correction term decreases. Thus, in the approximate solution 
we have the power 3(z ta — 0.6) term. Also we find that z ta ~ 0.6 in order to have z v - u ~ in the 
concordance model. 



Since we know the exact value of r c from Eq. (2.13), we are able to estimate the goodness of 
the approximate solution given in Eq. (2.24). The errors of r c for the different models obtained 
from Eq. (2.24) are shown in Fig. 4. We show the errors of r^ naly against with zt a = 1.1 for 
the different u>de models in the left panel of Fig. 4. The dotted, dashed, solid, and dot-dashed 
lines correspond to a; dc = —0.8, —0.9, —1.0, and —1.05, respectively. If we limit 0.24 < Sl^ < 0.3, 
then the errors are less than 1.6 % for all the models. We also depict the errors of r^ naly against 
z ta with = 0.27 for the different c^dc models in the right panel of Fig. 4. Again, the dotted, 
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Fi gure 4: Error (%) of the analytic solution of t c obtained from Eq. (2.24). a) Error of f^ 11 ^ 
against when z ta = l.l(z v i r ~ 0.3) for the different Wde = —0.8, —0.9, —1.0, and —1.05 models, 
denoted by the dotted, dashed, solid, and dot-dashed lines, respectively, b) Error of r* naly against 
Zta when = 0.27 for the different lu^ c models. 

dashed, solid, and dot-dashed lines correspond to u>de = —0.8, —0.9, —1.0, and —1.05, respectively. 
The errors of r analy for lu^c > — 1.0 are less than 1 % for the entire range of z^ a . 

We claim that Csk gives the better fit than £ ws for both general Wde and a wide range of the 
cosmological parameters. As we expect, both give a good approximate value when z ta increases. 
Again this comparison is possible because we know the exact collapsing time r c = 2rt a for each 
Wde model. In Fig. 5, we demonstrate the evolution of both numerical and approximate analytic 
solutions y for some cases. The circle and rectangular points represent the numerical y values when 
we use Csk and ( ws , respectively. The diamond point (the grey point at y = 0) indicates the exact 
value of t c obtained from Eq. (2.13). The dashed line depicts the approximate solution of y given 
in Eq. (2.24). In the left panel of Fig. 5, we show both numerical and analytic solutions when 
u>de = —0.9, fJ^ = 0.3, and z ta = 1-0. We also demonstrate them when Wd e = —1-0, = 0.3, 
and zta = 0.4 in the right panel of Fig. 5. Even though we show z ta = 0.4 case in this figure, we 
will show in the next section that there exists the minimum z ta for each model when we use the 
virial theorem. 



3 Virialization 

If only the matter virializes, then we are not able to use the energy conservation in the virial 
theorem because the dark energy does not virialize with the dark matter [8, 17, 18, 19]. 

^cluster ( z ta.) + t^dc( z ta) > ^ ^cluster ( z vir) + 2[/d e (£ v ir) > (3-1) 

where U c i uste r and U^e are the potential energies associated with the spherical mass overdensity 
and with the DE, respectively. Thus, in the presence of the smooth DE, the dark matter can 
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Figure 5: Comparing numerical solutions of y when we use Cws and Csk with y ana i y for the different 
models, a) r v.s. y when ujde = —0.9, = 0.3 and z ta = 1. b) r v.s. y when uj& e = —1.0, 
= 0.3, and 

z t& = 0.4. The circle and rectangular points represent the numerical y values when 
we use Csk and Cws, respectively. The diamond points at y = indicate the exact values of r c . 



reach a quasi-equilibrium state where virialization holds instantaneously. To estimate the virial 
radius y v i r = we use the assumption Pd e (z v h) = Pde^ta) as in the reference [17]. With this 
approximation, we are able to obtain the equation of y v - n = ^ L for any value of uj^c '■ 
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^ like that in the EdS universe 



An interesting feature in this equation is that we obtain y v \r 
when ui^e = — g for any given values of and Zt a [13]. Thus, we are able to obtain x v - 1T of the 
= — | model for the different values of and z ta from Eq. (2.18). With simple algebra by 
using Eqs. (2.13), (2.14), (2.18), and (2.19), we obtain 
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tain the analytic expression of x 
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in the 



this equation can be solved in a non-algebraic way. However, be represented by the 

approximate analytic form for the wide ranges of the cosmological parameters (fl^ and zta.) if one 
analyzes the form by interpolating x Yir values. 

For the EdS universe, Eqs. (2.13), (2.18), and (2.19) become 
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" ArcSinl^y] + \ = V((r ~ \) = f (x§ - l) , (3.5) 

C = ( 3 i?, (3-6) 

where we use F[^,— 1 — ^j-, 0] = 1 and r ta = §• The above equations can be generalized 
to any value of Wd c when approaches to 1 or z tSL goes to a very high value (for most of the 
parameter spaces z tSL > 5). In the EdS universe, y v ; r = ^ from Eq. (3.2). The commonly used 
assumption to obtain the nonlinear overdensity A c in the EdS universe is that r c = 2rt a = | 
which is the collapsing time for y c = even though one uses y v - ir = | to obtain A c . By using this 

2 

assumption, one obtain x c = 23 from Eq. (3.4). Thus, one obtains 



A »_ < ^ = - « =ft. aM . (3.7) 




However, we can obtain the correct value of A v ; r for the EdS universe because we know the exact 
relations between r, x, and y. After replacing y v - ir = \ into Eq. (3.5), we obtain the virialization 
time T v i r = + 1 < t c . From this correct value of r v ; r , we obtain x v \ T = + |)i by using Eq. 
(3.4). Thus, the correct value of the non-linear overdensity A v i r for the EdS universe is [13] 




As (Jde becomes large, the structure is formed earlier and A v ; r becomes large. There is only gravity 
involved in the structure formation in the EdS universe and thus it requires less overdensity than 
the one with the pressure. Thus, the minimum overdensity for general Wdc models is about 147 
instead of 178. We show this in Tables 1 and 2. As z ta increases, A v i r approaches to 147 instead 
of 178. 

The ratio of cluster to background density at the virialized epoch z v - ir for lo^c = — | becomes 
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^vir = ^v^virj 




18vr 



- X vir 



At p 13 5 X v j T \ 

*\ r 2' 2' 2' Qta > 



(3.9) 



where we use Eqs. (2.3), (2.6), and (2.19). We also use the fact that y v - 1T = | independent of Qta- 
We are able to obtain A v i r for the given cosmological parameters by using Eqs. (3.3) and (3.9). 
It is also useful to express z v - 1T by using zt a and x v ir from Eq. (2.9) 

z vir = i±^_l. (3.10) 

X vir 

Now we are able to extend the previous consideration to general values of uj& e by using Eqs. 
(2.13), (2.22), (2.24), and (3.2). We show several quantities for the different models in Tables 1 
and 2. We need to clarify several things in these tables. First, structures start to form earlier for 
lager oj^c and thus A v i r rises with increasing oj^c- Thus, we are able to distinguish between the 
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<^dc 


*ta = 0-7 


Zta = !•! 


^ = 0.24(z vir ~ 0.02) 


= 0.3(z vir ~ 0.04) 


= 0.24(z vir ~ 0.31) 


= 0.3(z vir ~ 0.33) 


Csk 




Uvir 


A v ir 


Csk 


•^vir 


2/vir 


A v ir 


Csk 




2/vir 


A v ir 


Csk 


•^vir 


UvlT 


A v ir 


1 

3 


11.2 


1.62 


0.50 


384 


9.8 


1.60 


0.50 


322 


10.2 


1.61 


0.50 


339 


9.0 


1.59 


0.50 


289 


-0.6 


9.1 


1.65 


0.53 


282 


8.2 


1.63 


0.52 


250 


8.0 


1.61 


0.52 


236 


7.3 


1.59 


0.52 


214 


-0.8 


8.1 


1.66 


0.53 


244 


7.4 


1.64 


0.53 


222 


7.1 


1.60 


0.52 


201 


6.7 


1.58 


0.52 


188 


-0.9 


7.8 


1.66 


0.54 


231 


7.2 


1.64 


0.53 


213 


6.8 


1.60 


0.53 


193 


6.5 


1.58 


0.52 


182 


-1.0 


7.6 


1.66 


0.54 


222 


7.1 


1.65 


0.53 


212 


6.6 


1.60 


0.53 


189 


6.4 


1.59 


0.52 


181 



Table 1: Csk) ^vin 2/vin an d A v ; r with the given values of the cosmological parameters for the 
different DE (t^de) models. As Wde becomes large, the structure is formed earlier and A v i r becomes 
large. 





Zba. = 1-4 


Zta. = 5.0 


n v m = 0.24(z vir ~ 0.53) 


W m = 0.3(z vir ~ 0.54) 


ny n = 0.24(z vir ~ 3.0) 


W m = 0.3(z v ir ~ 3.0) 


Csk 




Uvir 


A v ir 


Csk 




Uvir 


A v ir 


Csk 


X v ; r 


Vvir 


A v ir 


Csk 




2/vir 


A v ir 


1 

3 


9.6 


1.60 


0.50 


316 


8.6 


1.58 


0.50 


271 


7.2 


1.55 


0.50 


215 


6.8 


1.50 


0.50 


183 


-0.6 


7.4 


1.59 


0.52 


216 


6.9 


1.57 


0.51 


198 


5.9 


1.52 


0.50 


161 


5.8 


1.50 


0.50 


154 


-0.8 


6.7 


1.57 


0.52 


184 


6.4 


1.56 


0.52 


176 


5.7 


1.50 


0.50 


151 


5.6 


1.50 


0.50 


150 


-0.9 


6.4 


1.57 


0.52 


178 


6.2 


1.55 


0.51 


171 


5.6 


1.50 


0.50 


150 


5.6 


1.49 


0.50 


147 


-1.0 


6.3 


1.57 


0.52 


175 


6.1 


1.56 


0.51 


170 


5.6 


1.50 


0.50 


150 


5.6 


1.49 


0.50 


147 



Table 2: Csk> ^vir, Uvh, and A v ; r with the given values of the cosmological parameters for the dif- 
ferent DE (wdc) models. Independent of model parameters, A v ; r converges to 147 as z ta increases. 



models from the observed cluster at a given z v i r . The other important feature is that as and 
Zta increase, the model dependence of the observed quantity A v ; r decreases. Thus, the low-redshift 
clusters are the better objects to be used for studying the properties of DE. Again as zt a increases, 
Avir approaches to 147 instead of 187r 2 ~ 178. 



From Eq. (3.9), we are able to investigate the linear perturbation at early epoch A l+^iin 

A = i + Snl ™S 1 + <5 Iin = 1 + | (V^r) 1 = 1 + | ( v^^) § , (3.11) 
where we use H ta = in the early (matter-dominated) epoch. This is equal to the famous 

2 

result for the EdS universe (a/C = x an d ^ mta = 1), Su D = f(xt^) 3 [2, 13]. As r approaches to 
zero, the exact solutions of x and y given in Eqs. (2.13) and (2.17) become the good approximate 
solutions for general values of Wde- I n this limit, y v ; r — \ and thus T vir — (| + ^) T ta from 
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Eq. (2.17). Thus, 



) 



2 
3 



20 



+ 9vr)3 ~ 1.58 



(3.12) 



where we use a/C = The above result given in Eq. (3.12) is true for any value of lo^c because 
uJd c dependence of x and y is disappeared in the r — > limit. From this critical density threshold 
^iin(^vir), we are able to obtain 6\i n at any epoch by using the relation 



where D g is the linear growth factor. There is the exact analytic form of D g for the DE model 
with the constant eos [20]. 

From the analytic forms of dynamical quantities x, y, and (, we are able to estimate the 
abundances of both virialized and non-virialized clusters at any epoch. Also the temperature 
and luminosity functions are able to be computed at any epoch [21]. Thus, these analytic forms 
provide a very accurate and useful tool for probing the properties of DE. 

4 Conclusions 

The exact solution for the evolution of the background scale factor in a flat universe with the dark 
energy is given as a function of time. This provides the exact value of collapsing time and makes 
it possible to test the accuracies of both the analytic and the numerical solutions. We obtain the 
exact and approximate solutions of the evolution of the overdensity radius with high accuracy. 
From these solutions we are able to estimate the non-linear overdensity at any epoch. 

Even though the analytic forms of y and ( are obtained for the constant code models, they can 
be generalized to the slowly varying u>^ c > — 5 because they provide the very accurate values of 
them for the wide range of cosmological parameters. 

The virial theorem provides the exact value of the virial epoch for given models. From this, 
the value of the non-linear overdensity is obtained. In the high matter density and/or the high- 
redshift clusters, the model dependence of the observed quantities becomes weaker. We are able 
to categorize the physical quantities of the low and high redshift clusters for the given models. 

From the correct value of the virial radius of an EdS universe, we obtain the correct linear 
and non-linear overdensities for an EdS universe, which are lower than the known values. In this 
paper, we show that clusters are formed earlier than those from the conventional model. Thus, 
the mass of the cluster is larger than one obtained from the conventional model. This will affect to 
the prediction for the different models for the cluster abundances, the weak gravitational lensing, 
etc. 



DM 

Dg (^vir) 




(3.13) 
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A Appendix 



We derive the exact solution of r as a function of x given in Eq. (2.13). After replacing the 
variables Z = and T = (^)~ 3wdc , Eq. (2.13) becomes [16] 



Jo 



dx' 



_1 , x'-^Ao- 1 
^ Qta 



3 

3wd< 



r 1 _i l_ i 

/ t 2 "dc(i-zr)-2 

Jo 



2 3 

2dT = -x?F 



x 



-3u> de 



1 

2 ' 2tJdc ' 2Wde ' Qta. 

(A.l) 

where F is the hypergeometric function and we use the gamma function relation r[l + b] = bT[b]. 

The evolution equation of r as a function of y when lo^c = — | is given in Eq. (2.16). The 
exact analytic solution of this equation is known [13] : 



ArcSinfVy] - yy(l - y) = VCnT , when r < r ta i3 , (A.2) 
\/y(l - y) - ArcSin[^/y] + ^ = \/Ci3 {r ~ 7"tai3) , when r > r ta i3 , (A.3) 

where £13 is given in Eq. (2.19) and r ta i3 is T ta for uj^c = — 3- ftai3 is obtained from Eq. (2.14) as 



Ttal3 — -^F 



1 3 5 



1 



2'2'2' Q 



ta!3 



(A.4) 



where Q ta i3 is Qta for u Ae = -\. 
The equation governing y for Wdc 



■1-1 



-1 is given by Eq. (2.20) with a = ( + Q^a 



C(J"i) + cfe(y 2 -i) 



C 2 ± T . 



(A.5) 



The above equation (A.5) is solved analytically by [16] 
dy' 



r 



ta 



2(QtaC - 2) 



EllipticF ArcSin[B],C 



(3 + vi + 4g ta c) 



EllipticPi A,ArcSin[B],C 



where A 



2y/l + 4QtaC 



3 + VI + 4QtaC ' 



B = 



3 + 



(2QtaC~4) 



+ vi + 4g ta c 



V2(l + 4Q ta C) 



,c = 



2V1 + 4QtaC 



-i + 2Q ta c + vi + 4g ta c' 



(A.6) 



where EllipticF and EllipticPi represent the elliptic integral of the first kind and the incompelete 
elliptic integral of the third kind, respectively. 



We also check the linear perturbation at early epoch, A 



T->0 



1 + <5i in . We use Eq. (A.2) 



because we consider r —> 0. Even though this solution is exact only for lo^ c = — it is safe for 



3' 



us to use it for the general values of Wde because the second term in Eq. (2.15) becomes zero in 
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y — > limit. Thus, <5ii n obtained from this limit is independent of the value of uj<{ c . |x2 ~ r and 
|ya + Iy2 ~ y/Qr as r — > 0. Thus we obtain 

A=1 + ftn ._ < (£)'= !< ^( 1+ + ^ fe^-f^fr)'. (A.7) 

We thank S. Basilakos and D. F. Mota for useful comments. 
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